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1. Introduction: Massive (gauge-fixed) gluons and Gribov ambiguities 


Progress in the understanding of the phase structure of QCD is mainly driven by lattice results. 
It is, by now, well established that there is a crossover from a hadron phase to a quark-gluon 
plasma phase along the temperature axis |||, §]. In contrast, a simple perturbative description of 
the confinement-deconfinement transition seems out of reach even in the case of pure Yang-Mills 
theories because weak coupling techniques break down at temperatures of the order of the transition 
temperature. Continuum calculations are either based on nonperturbative approaches such as the 
functional renormalization group (FRG) or Dyson-Schwinger equations (DSE) ||3|, or on the 
use of phenomenological models [|>01- 

However, the failure of standard perturbative techniques at low temperatures should be inter¬ 
preted with care due to the Gribov problem. Indeed, the Faddeev-Popov (FP) quantization proce¬ 
dure ignores the Gribov ambiguities inherent to the issue of fixing a gauge in a continuous way in a 
nonabelian fheory [^] and is, af besf, a good approximation af high energies, where fhe Gribov am- 
biguifies play a negligible role. A fully consisfenf quantization procedure musf lake info accounl fhe 
Gribov problem. A well-esfablished example is fhe so-called minimal Fandau gauge in fhe confexl 
of lattice gauge-fixed calculations, which consisls in selecting one (arbifrary) Gribov copy on each 
gauge orbif. Fxplicif calculafions of fhe Euclidean vacuum propagators in fhis gauge have revealed 
fhaf fhe gluon behaves as a massive field in fhe regime of deep infrared momenla, whereas fhe 
ghosl remains massless |Q]. This demonslrales fhaf fhe BRST symmelry of fhe FP action—which 
prohibifs a gluon mass or, more precisely, a nonzero value of fhe inverse gluon propagator af zero 
momenfum—is nol realized in a fully gauge-fixed setting. In facl, if is also well-known fhaf lattice 
implemenfafions of fhe BRST symmelry fypically lead lo ill-defined zero over zero expressions for 


physical observables []10|]. 

Unforlunalely, fhe minimal Fandau gauge cannol be implemenled in terms of a local renor- 
malizable acfion and is nol suited for conlinuum approaches. An early attempl to go beyond fhe 
FP quantization consisls in reslricling fhe palh infegral over gauge fields lo fhe firsl Gribov region, 
where fhe FP operator is posifive definite sinii- The original proposal predicfs a vanishing gluon 


propagator af zero momenfum, which agrees wifh lattice resulls in fhe Coulomb gauge 012| ] bul 
nol wifh Ihose in fhe Fandau gauge menlioned above. A refined version, which includes various 
condensales of mass-dimension Iwo gives resulls in agreemenl wifh fhe lattice dala in fhe Fandau 
gauge [|T^]. Finite femperafure effecls have been invesligaled in fhis confexl in Refs. 10,111,11611. 

Yel anolher approach to fhe Gribov problem, advocaled in Refs. is lo acknowledge fhe 

facl fhaf fhe BRST symmelry is likely nol lo be realized in a consisfenf gauge fixing procedure and 
to consider fhe minimal extension of fhe FP acfion consisfenf wifh localily and renormalizabilily. In 
fhe case of fhe Fandau gauge, fhis is given by fhe Fandau limil of fhe Curci-Ferrari (CF) model [ [I^ 
which corresponds to fhe FP acfion augmented by a bare gluon mass term. One-loop calculafions 
of fhe vacuum ghosl and gluon propagators in Ibis model have been shown lo be in quanfilafive 
agreemenl wifh fhe aforemenfioned lattice dala [|l7|]. Moreover, Tissier and Wschebor have shown 
fhaf fhe model admits infrared-safe renormalization group trajectories, where the coupling constant 
remains finite at all scales, which allows one to push perturbative calculations down to deep infrared 
momenta. The gluon mass regulates the infrared divergences responsible for the Fandau pole in the 
FP theory. This approach has also been successfully applied to the one-loop calculation of three- 
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point vacuum correlators in Yang-Mills theories [20], to the vacuum propagators of QCD [21 ] and 
to the ghost and gluon propagators of Yang-Mills theories at finite temperature [^]. 

Finally, the relation of this phenomenological approach to the issue of Gribov ambiguities has 
been given a more solid theoretical foundation in Ref. [^^, where a novel quantization procedure 
for Yang-Mills theories was put forward, which consists in taking a particular average over all 
Gribov copies instead of trying to select a unique one. This avoids the Neuberger 0/0 problem 
mentioned above. For a suitably chosen averaging procedure, this can be formulated as a local ac¬ 
tion which is perturbatively renormalizable in four spacetime dimensions. In the case of the Landau 
gauge^ and for a suitable definition of the renormalized mass parameter, the resulting gauge-fixed 
theory has been shown to be perturbatively equivalent to the CF model for the calculation of ghost 
and gluon correlation functions. The bare gluon mass is related to the averaging weight which lifts 
the degeneracy between Gribov copies. 

The results mentioned above in the CF model show that this approach leads to a better-behaved 
perturbation theory than the usual FP procedure. A qualitative but instructive analogy can be made 
with the physics of the Kosterlitz-Thouless transition in the XY model in statistical physics [^]. In 
terms of usual spin waves, the physics of the transition is intrinsically nonperturbative: perturbation 
theory shows no sign of the transition at all orders. Alternatively, taking explicitly into account the 
vortex excitations of the XY model leads to a valid perturbative description of the transition. 

These considerations have motivated the works of Refs. [^, summarized below, where 
the confinement-deconfinement transition of pure SU(A) Yang-Mills theories at finite temperature 
is studied in the context of the modified (massive) perturbative approach described above. A crucial 
ingredient is to devise a loop expansion which correctly captures the Za? center symmetry of the 
theory, which is spontaneously broken in the deconfined phase. This is easily achieved in the 
context of background field methods. In particular, we consider the minimal massive extension 
of the standard Landau-DeWitt gauge, the background field generalization of the Landau gauge. 
A one-loop calculation of the effective potential for the Polyakov loop—the order parameter of 
the ZjM transition—correctly describes a confined phase at low temperatures and a deconfinement 
transition of second and of first order for the cases N = 2 and N = 3 respectively, with transition 
temperatures in qualitative agreement with known values [^]. Finally, we briefly review the work 
of Ref. [ 27 [ concerning the calculation of the critical temperature as well as of the thermodynamic 
pressure and entropy of the SU(2) theory at two-loop order. 


2. Massive Landau-DeWitt gauge 

We consider the Euclidean Yang-Mills action in d = 4 dimensions^ 

~ 2 J^ ) ( 2 . 1 ) 

where F^v = d^Ay — dyA^ — ig[A^,Ay], with g the coupling constant, and = A“t“, with SU(A) 
generators F normalized as tr(t"t^) = 5“^jl. Finally, = Jq dzf d^x, with p = IjT the inverse 
temperature. We quantize the theory using the background field method [ ^^ with a background 

* A more general class of (nonlinear) covariant gauges has been considered in Ref. [0]- 

^For simplicity, we ignore regularization and renormalization issues in this contribution. For details, see Ref. 
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field in the temporal direction and in the Cartan subalgebra spanned by the generators 

with k = . ,N —\. Our gauge-fixed action reads 


S = 





+ m^a^a^ 


+ ID^cD^c + lihD^a^ 


( 2 . 2 ) 


with = A^—Api the fluctuating gluon field, h a (real) Nakanishi-Lautrup field and c and c the 
FP ghost and antighost fields. Here, D^(p = d^(p — ig[A^,(p] and = d^(p — ig[Aj^,(p]. Apart 
from the bare gluon mass term oc m^, Eq. (^) is nothing but the FP action corresponding to the 
Landau-DeWitt gauge: = 0. The relation of the gluon mass term in Eq. (^) with the Gribov 

ambiguities of this gauge condition is discussed in Ref. [|2^. 

To evaluate physical observables at zero sources, one can minimize the following background 
field effective potential 

, rfA. m = ni 

(2.3) 




pa 

where = PgA^, Q. is the spatial volume, and where we have subtracted the zero temperature 
contribution Vyac- Here, r[A,(p] is the effective action in presence of the background, with cp = 
{a^,c,c,ih). The loop expansion of the potential (2.3) corresponds to an expansion in powers of g 
with gA() ~ ^(1). We write the corresponding series as 


E(r,r^) = ^EW(r,r^), (2.4) 

«>0 

with ~ the n-loop order contribution. The tree-level contribution vanishes identi¬ 

cally, V^^\T,r^) = 0, and the first nontrivial term is the one-loop contribution. 


3. One-loop results 

Remarkably, the leading-order, one-loop contribution (T, r^) already captures the essential 
physics of the transition. This arises from the interplay between the respective contributions of 
the massive (gluon) and the massless (ghost) degrees of freedom. The one-loop expression reads 

r') = ^ {iTr In - Tr In |, (3.1) 

with and A^^ the tree-level inverse propagators in the {a,h) and in the (c,c) sectors in the 
presence of the background field. The notation Tr on the right-hand side involves a trace over 
color and Eorentz indices as well as a sum over Matsubara frequencies and an integral over spatial 
momenta. The Matsubara sums are easily evaluated using standard contour integration techniques. 
Below we summarize the main results for the gauge groups SU(2) and SU(3) [^]. 

3.1 Second order phase transition for SU(2) 

In that case, the Cartan subalgebra has a single direction. Thanks to the symmetries of the 
problem one can reduce the analysis to the domain r G [0, n]. Introducing the function 

^ j r/^^^jln ^1 — -l-ln ^1 — 2e^^®''cosr^ (3.2) 
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Figure 1: The SU(2) dimensionless background field potential 'fiij^r) = V{T,r)lT^, normalized to its 
value at the confining point r — n, for temperatures T — Tc (black), T <Tc (blue), and T > Tc (red). The 
minimum is at r = tt for T <Tc (black line) and continuously moves toward zero for higher temperatures. 
The left figure shows the Weiss and inverted Weiss potential at high and low temperatures respectively. The 
right figure is a close-up view around Tc. 


which is such that, for r E [0,27r], 


^o(r,r) 


jt4 

~6 


[r-n) 

2n^ 


-{r-K)'^ + 


71 ^ 

To 


the one-loop background field potential can be written as 




(3.3) 


The first term on the right-hand side is the contribution from the massive gluons and the second one 
is due to the incomplete cancelation of the massless modes in the {a,h) and in the (c,c) sectors. At 
the same order of perturbation theory, the order parameter of the Z 2 transition, i.e., the average of 
the traced Polyakov loop is given by its tree-level expression 


£{T) = cos 




(3.4) 


where is the absolute minimum of the one-loop potential (^) at fixed femperafure T. In 

fhe high femperafure limif, fhe effecfive gluon mass is negligible and one recovers fhe well-known 


Weiss pofenfial |29], which corresponds fo fhe sfandard FP quanfizafion, 




(3.5) 


The minimum sifs af r = 0, corresponding fo a deconfined phase, where i{T) In confrasf, af 
low femperafures fhe confribufion from massive modes is exponentially suppressed and one obfains 
an inverfed Weiss pofenfial from fhe massless modes, as firsl discussed in Ref. [Q] : 

(3.6) 


The minimum is now at r = 71 , which corresponds fo a confined phase, wifh £{T) = 0. A de- 
failed sfudy of fhe expression (3.3) reveals a second order phase transition from the confined fo fhe 
deconfined phase af a critical femperafure Tc/m Ri 0.33, as illusfrafed in Fig. ^ 
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Figure 2: The SU(3) dimensionless background field potential 'P{T,r^,rs) =V{T,r^,r^)/T^ in the rg = 0 
direction (where the absolute minimum always sits), normalized to its value at the confining point (rg = 
47r/3,r8 = 0), for temperatures T — Tc (black), T < Tc (blue), and T > Tc (red). The right figure is a 
close-up view around Tc- One observes a finite jump of the absolute minimum atT = Tc. 



Figure 3: The Polyakov loop as a function of the temperature normalized to the transition temperature for 
SU(2) (blue) and SU(3) (red), from the leading-order calculation. 


To estimate we use the value of the effective gluon mass parameter obtained by fitting the 
lattice data for the vacuum gluon propagator in the Landau gauge against the tree-level expression 
in the present theory at vanishing background. This is motivated by the fact that, for T < Tc, 
the physical background is given by Aq = 7 lT jg and thus vanishes at T = 0. In J = 4, we obtain 
m = 710 MeV, which gives Tc ~ 238 MeV. The typical lattice value p0[ ] is = 295 MeV and the 
most recent value from FRG/DSE studies^ 0 is = 230 MeV. Although such a comparison is 
only qualitative because of the issue of properly setting the scale, it is remarkable that the present 
one-loop result falls in the right ballpark. 

3.2 First order phase transition for SU(3) 

The Cartan subalgebra has now two directions, conventionally taken as (rgjrg). The calcula¬ 
tion of the one-loop contribution ( |3.1D is, again, elementary. As in the previous case, it reduces 
to the perturbative SU(3) Weiss potential at high temperatures and to a confining, inverted Weiss 
pofenfial af low femperafures. The fransifion is now of firsf order, as shown in Fig. wifh a fran- 
sifion femperafure Tc/m k, 0.36. As before, we esfimafe fhe gluon mass parameter from fils of 
SU(3) laflice dafa for fhe vacuum propagafor in fhe Landau gauge. We find m = 510 MeV, which 

^The authors of Ref. [Q] mention that their FRG result is modified to = 300 MeV when some backreaction 
effects—neglected in their main study—are included. 
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Figure 4: Temperature dependence of the Polyakov loop at leading (red) and next-to-leading (blue) orders 
in the SU(2) theory. The horizontal dashed lines denote the corresponding asymptotic values at high tem¬ 
perature, denoted here by ilj and respectively. The respective one- and two-loop critical temperatures 
are indicated by vertical dashed lines. We have also indicated the temperature T* at which the leading-order 
Polyakov loop reaches its asymptotic value. 


yields = 185 MeV. With the same word of caution as in the SU(2) case, we mention the results 
from lattice calculations [^0|], = 270 MeV, and from FRG/DSE studies [Q], = 275 MeV. 

Our one-loop results are summarized in Fig. ^ which shows the temperature dependence of the 
Polyakov loop for both SU(2) and SU(3). We mention that all the continuum approaches men¬ 
tioned here produce a rapid increase of the Polyakov loop above Tc, in sharp contrast with the 
existing lattice data |1], 


4. Two-loop results 


An important aspect of the present perturbative approach is that the above leading-order re¬ 
sults can be systematically improved by computing higher order corrections. The next-to-leading- 
order, two-loop contribution V^^\T,r^) to the background effective potential has been computed in 
Ref. for the SU(2) theory. The calculation is considerably more involved than the leading-order 
one, but it can be done almost completely analytically, up to two-dimensional radial momentum 
integrals involving Bose-Einstein ditribution functions, which can be easily computed numerically. 

Such a calculation is of interest, first, in assessing the convergence of the present perturbative 
approach and, second, because the two-loop contributions resolve some unphysical artifacts of 
the leading-order results. The first aspect is illustrated in Fig. ^ which shows the temperature 
dependence of the Polyakov loop at one- and two-loop orders. For illustration, the next-to-leading- 
order expression reads 


i{T) =cos 


.{T) 


1 + 


g^m 


3271 


+ ^ T sm-" 




2/1 T 

where rmin(T) is the absolute minimum of the two-loop background field potential and 

' k^dk f 1 


t(T,r)= r 
Jo 




k^lel 


cosh (j3 A:) — COST cosh(j3eA:) — COST 


> 0 . 


(4.1) 


(4.2) 
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The critical temperature now depends on both the effective gluon mass and the coupling. We es¬ 
timate the latter using the same strategy as before: we fit the lattice data for the vacuum ghost 
and gluon propagators in the Landau gauge against the perturbative expressions at vanishing back¬ 
ground field at the appropriate order of approximation.^ We obtain a 20% correction for the critical 
temperature: Tc ~ 284 MeV. This is summarized in the following Table. 


Tc (MeV) 

one loop 

two loop 

lattice ||3^] 

FRG/DSE 0 

SU(2) 

238 

284 

295 

230 (300) 

SU(3) 

185 

work in progress 

270 

275 


As for the second aspect mentioned above, there are indeed various unphysical artifacts of 
the one-loop results which disappear at two loop. The first one is the fact that the leading-order 
Polyakov loop reaches its asymptotic high temperature value lH, = 1 at a finite temperature T*, 
see Fig. ^ which results in a spurious singularity in thermodynamic quantities, as discussed in 
Ref. [H]. Such a singularity is not present at two-loop order, where the Polyakov loop only 
reaches its high temperature value^ asymptotically. Another example concerns the calculation 
of thermodynamic quantities such as the pressure p{T) = —V{T,r^ia{T)) or the entropy density 
s{T) = dp/dJ ; see Fig. ^ At leading order, we find a narrow range of temperatures around 
where p and/or s are negative. Again, both unphysical behaviors are cured by the two-loop contri¬ 
butions which are important near the critical temperature, as one might expect. This is discussed in 
detail in Ref. |27]. 

It is also interesting to compare the present calculation of the pressure with a nontrivial back¬ 
ground field to that at vanishing background, which is equivalent to the (massive) Landau gauge. 
Figure ^ shows that in the latter case, the leading-order pressure is always negative below a cer¬ 
tain temperature^ and that the two-loop corrections make the problem even worth since the pressure 
seems always negative. Clearly, neglecting the nontrivial background in the low temperature regime 
amounts to expand around an unstable point. This suggests that the strong systematic effects ob¬ 
served in lattice calculations of the Landau gauge gluon propagator near criticality [ M ] might be 
resolved by performing lattice calculations in the Landau-DeWitt gauge, e.g., along the lines of 
Ref. 

Finally, we mention that there remain unphysical contributions to thermodynamic quantities 
from the massless (e.g., ghost) modes of the gauge-fixed theory. For instance, these leads to spu¬ 
rious contributions to the pressure at low temperatures, as illustrated in Fig. If, as we have 
discussed above, these modes play a crucial role in ensuring a confining potential for the Polyakov 
loop at low temperatures, they should not directly contribute to physical quantities such as the 
pressure or the entropy. At high temperatures, the contributions from unphysical modes cancel out 
thanks to the effective restoration of the BRST symmetry. However, the cancelation is only incom¬ 
plete at low temperatures, where the effective gluon mass term cannot be neglected. This important 


^We employ the same renormalization conditions as in Ref. |p^, where such fits have been performed. The best fits 
are obtained with m = 680 MeV and g = 1.5 at a scale ft = 1 GeV. As a rough error estimation, we have checked that a 
30% change in the value of the coupling results in a 10% change in the critical temperature. 

^The high temperature limit considered here is only formal in the sense that it ignores important physical effects 
such as the running of the coupling or the physics of hard thermal loops. 

similar observation has been made in a different context in Ref. |p^. 
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p (GeV-'*) 



p/T^ 



Figure 5: Left; Thermodynamic pressure at one-loop (red) and two-loop (blue) orders, obtained from the 
minimum of the background field potential as a function of the temperature. The plot in the inset is a zoom 
on the low temperature region. The respective one- and two-loop critical temperatures are indicated by 
vertical dashed lines. Right; Comparison between one-loop (red) and two-loop (blue) results in the massive 
Landau-DeWitt gauge (plain) and in the massive Landau gauge (dashed). 


open conceptual issue requires further investigation. An interesting proposal in this direction has 
been made in Ref. [^] in the context of the Gribov-Zwanziger approach. 


5. Conclusions and perspectives 

In conclusion, we have proposed a modified perturbative approach to the physics of the static 
quark confinement-deconfinement transition that takes into account the Gribov issue in (gauge- 
fixed) continuum calculations in an effective way. This is based on a massive extension of the 
Landau-DeWitt gauge in the context of background field methods. The approach correctly captures 
the physics of the phase transition in Yang-Mills theories already at leading order in perturbation 
theory and next-to-leading-order corrections clearly improve the results. Although there are still 
important open questions, in particular, concerning thermodynamic quantities at low temperatures, 
we believe this is an interesting step towards a (semi)analytical understanding of the dynamics at 
work in the transition region. 

There are obvious extensions of the work presented here, among which, the calculation of the 
two-loop contributions in the SU(3) theory, the application of our approach at finite (real or imag¬ 
inary) chemical potential, or the calculation of the finite temperature ghost and gluon propagators 
at one-loop order in presence of the nontrivial background field. These are works in progress. 
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